clear; clc;

%% Generate the Baseline Results
load('Result_Main.mat');
StdVec      =   [1 1 1 -1 -1 60 -10]'*0.01/4;
                     
PP          =   Setup_PP(struct('KAPPA',0.00/4,'XI',1/82.5/4,'Pir_SS',0.02/4,...
                                'Flag_dE',0,'Flag_Price',0,...
                                'TrProb_H',0.99,'Prob_H',0.37,...
                                'TrProb_ext',0.99,'Prob_ext',0.33,...
                                'Fin_DomShare',0.79,'Fin_AdjCost',0.8,...
                                'Cons_FracT',0.3253,'Cons_FracH',0.60,...
                                'B_Total_SS',0.95,'b_lb',-0.20,...
                                'O_AdjCostVec_FI',[5.43;4.47],'O_AdjCostVec_RI',[5.24;4.51],...
                                'RHO_M_dom',0.68,'RHO_M_ext',0.81,'RHO_Y_H',0.50,'RHO_p_F',0.90,...
                                'Taylor_Pi',1.1,'Taylor_ir',0.87,...
                                'Cons_ElasTN',6.19,'Cons_ElasHF',6.19));

%% Different Results
TempPP              =   PP;
TempPP.Fin_AdjCost  =   0.8;
TempPP.Cons_ElasHF  =   6.19;
TempPP.Taylor_ir    =   0.95;

[IRF_New,SOLUTION_New] ...
                =   TempFun_NewIRF(TempPP,SS,MODEL,EquJac,StdVec);
            
Mom_Model       =   SubFun_IRF2Mom(PP,SS,MODEL,IRF_New);

Figure_Unit('TableGraphs//TempCalibration//',SS,IRF_New);
